Topological phases of dipolar particles in elongated Wannier orbitals 
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We show that topological phases with fractional excitations can occur in two-dimensional ultracold 
dipolar gases on a particular class of optical lattices. Due to the dipolar interaction and lattice 
confinement, a quantum dimer model emerges naturally as the effective theory describing the low- 
energy behaviors of these systems under well-controlled approximations. The desired hierarchy of 
interaction energy scales is achieved by controlling the anisotropy of the orbital dimers and the 
dipole moments of particles. Experimental realization and detection of various phases are discussed, 
as well as the possible relevance for quantum computation. 
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Unlike the ordered phase described by the Landau the- 
ory of phase transitions, a topological state is character- 
ized by its response to the changing of topology of the 
underlying manifold, instead of a local order parameter 
associated with symmetry breaking [l[. Theoretically, 
these phases are described by topological held theories, in 
which the low-energy elementary excitations often carry 
fractional quantum numbers of the constituent particles, 
a phenomenon known as fractionalization. These frac- 
tional excitations offer a concrete setting from which new 
physics emerges, e.g^, the realization of topological quan- 
tum computation ( 2] and reference therein). 

The most celebrated topological phases are perhaps 
the fractional quantum Hall (FQH) states described by 
the topological Chern-Simons gauge theory, whose low- 
energy excitations carry fractional charges and statis- 
tics of an electron. Another important example is the 
Zi topological phase in two-dimensional (2D) quantum 
dimer models (QDM) of frustrated spins, introduced by 
Rokhsar and Kivelson (RK) This novel phase is 

known as the resonating-valence-bond (RVB) state de- 
scribed by a topological Ising gauge theory and gives rise 
to deconfined fractional excitations such as spinons which 
carry spin 1 /2 but zero charge 

In ultracold atomic/molecular systems, efforts have 
been made in realizing various topological phases. For ex- 
ample, states similar to the FQH states were proposed in 
rotating atomic gases 111. There also have been sugges- 
tions to simulate QDM [12J and the RVB phase with ul- 
tracold atomic systems, but with the challenge of achiev- 
ing the special nonlocal interaction as pointed out 12 1. 
In fact, in contrast to its theoretical appeal, similar diffi- 
culties appear in most attempts to realize QDM and the 
RVB phase, for example, in the frustrated- Ising model 
[T3L [lj] and the easy-axis antiferromagnet [HI, [ly] • 

In this letter, we propose a different approach to re- 
alize and study the topological phases using ultracold 



dipolar particles (atoms or molecules) in specifically de- 
signed, non-standard optical lattices, in part motivated 
by the recent experiments of polar molecules ([lTHlij] and 
references therein). We exploit two novel features: 1) 
the orbital degree of freedom controlled by optical lat- 
tice confinement, and 2) the long-range interaction asso- 
ciated with dipoles. Using a site-to-bond duality map- 
ping [20l |. we show that the orbital motions on some op- 
tical lattices with dipolar interactions can be mapped 
onto a QDM: a particle in the elongated Wannier or- 
bital on the optical lattice corresponds to a dimer on the 
dual lattice; and the dipolar interaction between parti- 
cles enforces the "hard-core constraint" of the dimers. By 
tuning the optical lattice and the external field (the lat- 
ter aligns the dipole moments), this system shows great 
tunablity in exploring different phases of the QDM. For 
bosonic particles, within experimentally reasonable pa- 
rameter ranges, we shall show that various phases can be 
stabilized, including conventional crystalline phases and 
a Zi topological phase. The latter is similar to the RVB 
phase of the QDM, and will be referred to as the orbital- 
RVB (ORVB). Experimental detection of these phases 
and the fractional excitations in the ORVB phase shall be 
discussed. More interestingly, the orbital dimers in our 
model system can be either bosonic or fermionic, depend- 
ing on the statistics of the constituent particles, whereas 
the frustrated spin system only allows bosonic dimers. 
This opens a new possibility for studying fermionic QDM 
and even the mixture of dimers with different statistics. 

Orbital dimers in 2D optical lattices. We focus on 
a concrete example, a 2D Kagome lattice constituted 
by anisotropic potential wells [Fig. []Jb)]. Due to this 
anisotropy, the Wannier wavefunction at each site would 
be elongated and would be considered as a "dimer" as 
shown below. Note that arbitrary lattice potential can be 
generated in 2D for example by projecting a holographic 
mask through an imaging system [21] . Such potential 
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FIG. 1: (Color online) a) The Kagome lattice (solid lines) and 
its dual triangular lattice (dash lines). The green ellipse rep- 
resents the elongated wavefunction at each site of the Kagome 
lattice (the bond of the triangular lattice), b) Contour plot of 
the 2D double- well optical lattice potential (with I2 = 27i). 
Dipolar atoms (or molecules) reside at the intensity minima, 
which form the Kagome lattice, c) The local double-well po- 
tential and the wavefunction, which can be approximated by 
a Guassian. 



can also be achieved by independent sets of laser beams 
to produce an intensity pattern of (similar to Ref. [22I ]) 
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I 2 J2 cos 2 ( qj -r), (1) 

j=l,2,3 



where the amplitude of the wavevectors projected on the 
plane satisfy |qj| = v3|kj| (the z component is irrel- 
evant in a 2D system) and the azimuth angle of k^s 
(cfcs) are ki = in/ '3 (6^ = (j + l/2)vr/3). This in- 
tensity pattern induces an optical lattice consisting of 
double- well- like potentials [Fig. He)]. However, since 
the energy barrier between the two potential wells at one 
site is small, the double-well structure can be treated as 
a single anisotropic potential well, unless the energy of 
the trapped state is much smaller than the recoil energy, 
which is the limit we would avoid in this letter. 

In fact, for our interest, the precise shape of the po- 
tential is of little importance, as long as the ground state 
wavefunction at each site is sufficiently elongated and 
pointing towards the center of each hexagon. We call 
this elongated wavefunction an "orbital dimer" , because 
its spatial symmetry mimics that of a link variable de- 
fined on bonds. To make this analogy precise, we per- 
form a site-to-bond duality transformation. As shown in 
Fig. [TJa), under this mapping, an atom/molecule local- 
ized at a site of the original Kagome lattice (solid lines) 
becomes a dimer residing on a bond of the dual triangular 
lattice (dashed lines). Accordingly, particles on the orig- 
inal lattice are mapped to dimers on the dual lattice [23 1 . 

Dipolar interactions. If the dipolar particles are po- 
larized in the z-direction by an external field, their in- 
teraction is given by the long-range potential Vdd(r) = 
d 2 (l — 3z 2 /r 2 )/r 3 . Here, d is the dipole moment, r is 
the distance between the two dipoles, and z is the z- 



component of r. For a quasi-2£> system, the dipolar in- 
teraction is repulsive near the xy-plane, which leads to 
energy penalties for particles close to each other. The 
largest interaction is the on-site repulsion U and the 
next one is the repulsion between particles on the same 
hexagon of the original lattice, or cquivalently dimers 
joining at the same site on the dual lattice. Three such 
configurations i = 1, 2, 3 are shown in Fig.[2Ia)-(c). Even 
smaller is the interaction for parallel ( Vj|) a nd staggered 
(Vs) dimers, shown in Figs. 2(d) and 2(e) respectively. 
The effect of longer-range interactions will be studied be- 
low. In the limit U > Vi ^> V\\ > Vs, below a critical fill- 
ing (1/6 for the Kagome lattice), the large energy cost Vi 
(i = 1, 2, 3) forbids dimers to share a common site on the 
dual lattice. This enforces the "hard-core constraint" , a 
key feature of QDMs. At the critical filling, the dual lat- 
tice is "fully covered" by dimers. Namely, among the six 
bonds connected to any site, one and only one bond is 
occupied by a dimer. 

We now show that the desired hierarchy of interac- 
tions given above can be satisfied by controlling the dipo- 
lar interactions and increasing the elongation of Wannier 
wavefunctions. For this purpose, we approximate the 
Wannier function at each site of the Kagome lattice by 
these Gaussian "atomic orbitals" , 



w(x, y, z) = (af3j) 1/2 7r 3/4 exp[ 



x y g_ 1 
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with a > (5 ^> 7. Here, a strong laser field in the z 
direction is used to confine the particles within the xy 
plane (7 being small) and the x axis is chosen to be along 
the direction of the double- well potential as shown in Fig. 
[IJc). This simple ansatz would not provide an accurate 
value for U and Vs, but it suffices to estimate their order 
of magnitude and demonstrate the proper limit in which 
our system could be mapped to a QDM. 
Within this ansatz, we find 



U 



E{l-a 2 /t3 2 ) 



3a/?7 



a 2 /3 



where E(x) is the complete elliptic integral of the second 
kind. For tight confinement (small (3 and 7), U is exceed- 
ingly large, so it disfavors multi-occupancy of the same 
lattice site. On the dual lattice, this means that there can 
be at most one dimer on each bond. It is also straightfor- 
ward to compute the corresponding interaction Vi. For 
example, 



r 3 V 7T 



ddsec 3 9e 2 (y 2 - y 2 ) + 
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where y{&) = r 2 /(a 2 + /3 2 tan 2 (9), and r is the distance 
between the center of two dimers. 

Since the lattice structure is unspecified and can vary 
arbitrarily [2l|, we can treat a, j3, and 7 as free control- 
ling parameters and show that the limit U > Vi ^> V» > 
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FIG. 2: (Color online) Interactions between dipolar particles 
on the dual triangular lattice shown in Fig. [1] 

Vs can be achieved with elongated Wannier functions. 
First, it is useful to consider the case of isotropic Wannier 
functions, a = (3, and small wavefunction overlap, a <C r 
with r being the distance between two lattice sites. Then, 
the interaction energy is reduced to d 2 /r 3 , and the ratio 
Vi : V 2 : V 3 : Vj| : V s is 8 : 1.54 : 1 : 1 : 0.43. In particular, 
we see that Vz and Vj| are of the same order. A key point 
of our proposal is that by using elongated Wannier func- 
tions, the dipolar interaction for configuration i = 2, 3 
can be boosted with respect to the parallel configura- 
tion, since they have much larger wave function overlap 
(schematically shown in Fig. [5]). For example, for a = 0.3 
and (3 = 0.1 (in the units of the lattice constant), the ra- 
tio Vi : V 2 : V 5 : V\\ : V s becomes 30 : 8.6 : 8.6 : 1 : 0.49 
for 7 = /3/10 and 21 : 5.8 : 5.6 : 1 : 0.49 for 7 = ,8/5. 
Here V3 is almost one order of magnitude larger than 
Vj| and their ratio can be further enhanced by increasing 
the elongation of the Wannier function (through adjust- 
ing the optical lattice) and/or reducing 7. Notice that 
7 can be tuned independently. For a lattice with suffi- 
cient elongated Wannier functions, the mapping to QMD 
becomes asymptotically accurate as 7 is reduced. 

The orbital dimer model. At the critical filling, the low- 
energy fluctuations of the system are described by the 
following Hamiltonian after projecting out high-energy 
degrees of freedom (at the scale of U and Vi), 

H = -t eS J2(\ ~-)(W\+h.c.) 

+ Ve ff ^(| -)( -| + |W)(W|)+Vi fl , (2) 

where the summation is over all flippable plaquettes [5[ 
and Vlr describes higher order terms due to longer-range 
interactions (beyond those shown in Fig. [5]). Here the 
dimer has the same statistical property as the constituent 
particle, so it could be either bosonic or fermionic. At 
the critical filling, the dual triangular lattice is fully cov- 
ered by the orbital dimers. All local density fluctua- 
tions are gapped out due to large energy penalty of Vi 
(i = 1,2,3). Hence density remains homogeneous on 
large length scales and the dominant effect of longer- 
range repulsive interaction is to contribute a constant 
shift in the free energy. As a leading order approxima- 
tion, we will drop Vlr in the low-energy effective Hamil- 
tonian Eq. ^ and expect the orbital dimer model cap- 
tures the essential low-energy physics of the system. 




FIG. 3: (Color online) The phase diagram for the bosonic par- 
ticles and the particle configuration in the crystalline phases. 



At the critical filling, particle hopping gives rise to 
overlapping dimers (sharing the same site), which costs 
energy of order Vi . Such virtual hopping processes lead 
to quantum resonance of orbital dimers around a pla- 
quette (t e s), as well as corrections to Vj| and Vs through 
second order perturbations, of order t 2 /Vi. To the lead- 
ing approximation, we use one (average) energy scale t 
to estimate the amplitudes of relevant hoppings. The 
repulsive interaction between two parallel dimers, V c g, 
is set by the energy scale V\\ — Vs up to corrections of 
0(t 2 /Vi). By varying t and the interacting strength, t e g 
and V c s can be tuned to explore various parameter ranges 
of the QDM, except for the limit of V c g <C — t 2 /Vi, since 
V\\ > Vs for dipolar interaction. To estimate the rele- 
vant energy scales, we take d = 0.1 Debye, a = 0.25/xm, 
Vi ~ 4Vi>, 3 ~ 20Vj|, and t ~ lOnK. Then, V cS ~ V\\-V s ~ 
2.3nK and t cS - 2i 2 /Vi ~ 2.16nK. This shows that the 
condition U, Vi ^> V c s ~ t c g ^> T can be satisfied. 

For bosonic particles, the phase diagram at the crit- 
ical filling is shown in Fig. [31 using the knowledge of 
QDM [5]. As V cS /t ff increases, four phases appear con- 
sequently: the columnar phase, the y/l2 x \/l2 plaquc- 
tte phase (V e s/t e g ~ 0), the ORVB topological phase 
(Veff/teff < 1), and the staggered phase (V eS /t cS > 1). 
In the staggered and the columnar phases, the dimers 
are frozen in the ground state, while in the VT2 x Vl2 
state, they can resonate in each plaquette marked by the 
thick lines in Fig. [31 In terms of dipolar particles, all 
these three phases correspond to various commensurate 
charge-density-wave states, which break the point group 
symmetry of the underlying lattice. On the other hand, 
all the spatial symmetries of the lattice are preserved in 
the ORVB phase. However, this state has critical depen- 
dence on the topology of the underlying manifold. On 
a topologically nontrivial manifold, it shows degenerate 
ground states, topologically distinct from one another, 
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known as topological degeneracy. As a Z 2 topological 
phase, the ORVB phase has 4 degenerate ground states 
distinguished by two Z 2 topological indices if placed on a 
torus. As discussed below, the phenomenon of topolog- 
ical degeneracy not only provides unique experimental 
signatures for the ORVB phase, but also leads to the 
possibility of realizing quantum memory. 

In the context of frustrated spins the RVB state has 
two different types of fractional excitations: spinons (car- 
rying spin 1/2) and holons (carrying charge e), known as 
the spin-charge separation. For an ORVB state, how- 
ever, the orbital dimers are not spin singlets, and hence 
spinon excitations are forbidden. As a result, only one 
type of fractional excitations, holons, emerges. A holon 
carries only half the quantum number of an underlying 
dipolar particle, e.g., a half dipole, while the statistics of 
a holon may depend on microscopic details of the system 
26]. This phenomenon shares strong similarities with 
the deconfined monopole excitations in spin ice [27j . 

Using different optical lattices of double-well poten- 
tials, QDM on other lattices can also be achieved through 
similar constructions. For example, the following laser 
intensity results in a square-lattice QDM, 



! (Vr)+/ 2 



cos 2 (q.,- 



r). (3) 



with the in-plane component satisfying |qj| = 3|kj|/-\/2 
and the azimuth angles: = (2i + 1)tt/4 and # qj = 
2j7r/4. In this case, the ORVB phase shrinks into a RK 
point as expected for all bipartite lattices, while the crys- 
talline phases remain to occur in regions with dimer con- 
figurations slightly different from the triangular lattice. 

Experimental detection. Each of the phases discussed 
above displays a distinct interference pattern in Bragg 
scattering experiments. In the lattice shown in Fig. [T] 
for example, the ORVB phase has Bragg peaks at 
wavevectors commensurate with the reciprocal vectors 
of the lattice: (2tt, 2tt/v / 3) and (0,An/V3). In the stag- 
gered phase, however, the Bragg peaks are located at 
(mir, 2n7r/v / 3) for any integer m and n. Similarly, the 
columnar phase [Bragg peaks commensurate with vectors 
(tt, v^tt) and (0, 4tt/v^)] and the Vl2 x y/\2 phase [with 
vectors (tt/3, ±7r/\/3)] can also be distinguished. The 
six-fold rotational symmetry is broken in the staggered 
and columnar phase. This can be used to distinguished 
them from other phases in time-of-flight experiments. 

As in the experimental search for spin liquids, lack of 
symmetry breaking at low temperatures is an evidence 
of the ORVB phase in our system. One way to distin- 
guish it from the topologically trivial state is to detect 
the topological degeneracy in, for example, the "annulus" 
geometry, which can be achieved by punching a hole on 
the 2D optical lattice using an extra laser beam to form 
a forbidden region with high energy barrier. For a large 
enough hole, the ORVB ground state is doubly degen- 
erate. The two topologically-distinct degenerate states 



can serve as a "qubit" [28|. It can be "flipped" by cre- 
ating a pair of holons and then annihilating them after 
moving one holon around the hole of the annulus. Fur- 
thermore, topological quantum memory can be realized 
by punching multiple holes and manipulating the holons. 
For a hole of finite size, tunneling lifts the degeneracy and 
results in two nearly-degenerate states. The energy split- 
ting increases upon reducing the size of the hole. The size 
dependence of the energy splitting can serve as a direct 
evidence of Z2 topological order. However, it is challeng- 
ing to measure the splitting, which is on the order of pK , 
experimentally. Another distinct feature of the ORVB 
state is the holons. How to efficiently probe such frac- 
tional excitations remains an important open problem. 
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